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Abstract. I review the Globular Cluster Luminosity Function (GCLF) with emphasis 
on recent observational data and theoretical progress. As is well known, the turn-over 
magnitude (TOM) is a good distance indicator for early-type galaxies within the limits 
set by data quality and sufficient number of objects. A comparison with distances 
derived from surface brightness fluctuations with the available TOMs in the V-band 
reveals, however, many discrepant cases. These cases often violate the condition that 
the TOM should only be used as a distance indicator in old globular cluster systems. 
The existence of intermediate age-populations in early-type galaxies likely is the cause 
of many of these discrepancies. The connection between the luminosity functions of 
young and old cluster systems is discussed on the basis of modelling the dynamical 
evolution of cluster systems. Finally, I briefly present the current ideas of why such a 
universal structure as the GCLF exists. 



1 Introduction: What is the Globular Cluster Luminosity 
Function? 

Since the era of Shaplcy, who first explored the size of the Galaxy, the distances 
to globular clusters often set landmarks in establishing first the galactic, then 
the extragalactic distance scale. Among the methods which have been developed 
to determine the distances of early-type galaxies, the usage of globular clusters 
is one of the oldest, if not the oldest. Baum [5] first compared the brightness of 
the brightest globular clusters in M87 to those of M31. With the observational 
technology improving it became possible to reach fainter globular clusters and 
soon the conjecture was raised that the distribution of absolute magnitudes of 
globular clusters in a globular cluster system exhibits a universal shape, which 
can be well approximated by a Gaussian: 
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where dN is the number of globular clusters in an apparent magnitude bin dm, 
Too is the "Turn-Over Magnitude" and a the width of the Gaussian distribution. 
Also a representation by a "t5-function" 
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introduced by Seeker 82 found a wide-spread application. 

This distribution is called the "Globular Cluster Luminosity Function". In 
the following "Globular Cluster Luminosity Function" is abbreviated by GCLF, 
"Turn-Over Magnitude" by TOM, and "Globular Cluster System" by GCS. 

The conjecture that the GCLF is very similar in different galaxies, in par- 
ticular that the absolute magnitude of the TOM has an almost universal value, 
has been first suggested by Hanes [2Zj (but also see the references in this pa- 
per). The reviews of Harris & Racine [EH] , Hanes [EE|, Harris 0Ojj Jacoby et al. 
|3r>] . Ashman & Zepf 0, Whitmore Tammann & Sandage [EE], and Harris 

demonstrated both the solidity and the limitations of this conjecture. They 
also show the progress which has been achieved during the past 20 years both 
in terms of the number of investigated GCSs and the accuracy of an absolute 
calibration. 

The GCLF as a distance indicator has seen little application to spiral galaxies 
for several reasons: their GCSs are distinctly poorer than those of giant ellipti- 
cals, the identifications of clusters is rendered more difficult by the projection 
onto the disk, and the presence of dust causes inhomogeneous extinction. The 
investigation of GCSs of ellipticals or SO-galaxies is much easier due to the homo- 
geneous light background, the richness, and the absence of internal extinction. 

The application of GCLFs as distance indicators for early-type galaxies has a 
simple recipe: Given appropriately deep photometry of the host galaxy, identify 
GC candidates, as many as you can. In most cases this has to be done statistically 
by considering only objects with GC-like colors and by subtracting a hopefully 
well determined background of sources. Then measure their apparent magni- 
tudes, draw a histogram and fit a suitable function, for instance a Gaussian, to 
determine the apparent TOM. Under the assumption that every GCS has the 
same absolute TOM, one can use the galactic system and/or the M31 system 
to calibrate it in terms of absolute magnitudes. Real data, however, make the 
derivation of the TOM somewhat more difficult, which will be discussed below. 
The most important restriction is probably that we can follow the GCLF down 
to faint clusters only in two galaxies, the Milky Way and Andromeda. 

There is no physical reason why the GCLF should be a Gaussian or a t5- 
function. On the contrary, we shall later on learn about physical reasons why it is 
not a Gaussian. Closer scrutiny of the galactic cluster system indeed shows that 
its GCLF is not symmetric in that it exhibits an extended wing beyond the TOM 
for smaller masses (for example see Fig. 2 of Fall & Zhang ^H])- But Gaussians 
empirically are fair descriptions for the bright side and most observations of 
GCLFs of distant galaxies seldomly reach more than 1 mag beyond their TOMs, 
so this asymmetry is not relevant for measuring the TOM. 

Fig. shows the GCLF for the galactic system (upper panel) . It has been 
constructed on the basis of the " McMaster-catalog" (Harris [321 ) using the hori- 
zontal branch (HB) brightness as the distance indicator and adopting the relation 
Mv(HB) = 0.2 • [Fe/H] + 0.89 (Demarque et al. 10). The Gaussian fit results 
in M v = -7.56 ± 0.12 for the TOM and 1.2 ± 0.1 for the dispersion of the 
Gaussian. The lower panel shows a histogram of the masses, assuming M/L = 
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Fig. 1. The upper panel shows the globular cluster luminosity function for the galactic 
system together with a fitted Gaussian to the distribution. Only clusters with redden- 
ing E(B-V) less than 0.8 mag, and absolute magnitudes brighter than -4 have been 
considered in the fit. The lower panel shows the mass distribution in linear mass bins. 
The mass corresponding to the TOM is indicated. 



2, where the mass which corresponds to the TOM is indicated. In this linearly 
binned histogram, there is no striking feature at this mass. Indeed, the existence 
of a TOM is a consequence of the logarithmic magnitude scale in combination 
with a change of the power-law slope of the mass function. We come back to this 
in a later section. 

Throughout this review, we shall consider only TOMs in the V-band, be- 
cause most modern published data, particularly those from the Hubble Space 
Telescope, have been obtained in V. Other photometric systems, most notably 
the Washington photometric system (Geisler et al. [2E], Ostrov et al. [72], Dirsch 
et al. ^31) have been used for the investigations of GCSs as well. 

Given the previous excellent reviews on GCLFs as distance indicators, what 
can be the scope of this contribution? A lot of new data has been published 
during the last years and it is now possible to compare GCLFs of early-type 
galaxies with other distance indicators on the basis of a much larger sample 



4 Tom Richtler 



than has been possible before. The outstanding publication here is the catalog 
of distances based on surface brightness fluctuations (SBFs) (Tonry |H2)- We 
shall see that the absolute calibration of GCLFs indeed agrees very well with 
that of SBFs, demonstrating that most GCSs of elliptical or SO-galaxies show 
absolute TOMs which arc not distinguishable within the uncertainties of the 
measurements. Nevertheless, many discrepancies between GCLF distances and 
SBF distances exist. These galaxies are of particular interest and we shall discuss 
them as well. 

Beyond the usefulness of the GCLF as a distance indicator is the question 
why there exists such a remarkably universal structure. Is there a universal 
formation law for globular clusters, which operates in the same way in such 
different galaxies as the Milky Way and giant ellipticals? This problem has to 
do with the initial mass function of globular clusters and the evolution of GCSs. 
Much progress has been achieved during the last years, on which we will also 
report. 

2 Sources of uncertainty 

To begin with difficulties: Even if we trust the universal TOM, the actual mea- 
surement may appear straightforward according to the above recipe, but never- 
theless one encounters many sources of uncertainty. The identification of GCs 
as resolved objects from the ground is only possible for the nearest early-type 
galaxies, for example NGC 5128 (Rejkuba |77|1. Unfortunately, no modern in- 
vestigation of the GCLF of NGC 5128 exists until now. Observations with the 
Hubble Space Telescope can resolve the largest clusters in galaxies as distant as 
about 20 Mpc (e.g. Kundu & Whitmore HO], EH- Larsen et al. (23)- Therefore, 
the identification by ground-based observations normally has to use color criteria 
and the statistical discrimination against a " background" , which actually may 
consist of foreground stars and of unresolved background galaxies. That the lat- 
ter contamination is a strong function of the color system used, is nicely shown 
in Fig- El taken from Dirsch et al. ^Sjj which compares the color magnitude dia- 
grams V-I and Washington C-Tl for the GCS of NGC 1399, the central galaxy in 
the Fornax cluster. The C-Tl color recognizes many background galaxies, which 
have an excess flux in the blue band C, while they are not noticable in V-I, a 
color which has been widely used in HST investigations. 

The nearest large galaxy clusters (Virgo and Fornax) have distance moduli 
of about 31. Thus, TOMs for GCSs at this distance and beyond will generally be 
fainter than V=23.5, where the photometric incompleteness (depending on the 
data quality) plays an increasingly dominant role. Then there are the factors of 
the numbers of found GCs and the distance itself: If the photometry does not 
reach the TOM and/or the GCS is not very rich, the resulting TOM is naturally 
less well defined than in the case of a nearby, rich GCS. Then it may depend on 
the adopted shape of the luminosity function. However, empirically the derived 
TOM is not very sensitive to whether a t5-function or a Gaussian is used, at 
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Fig. 2. This plot has been taken from Dirsch et al. |13) . It shows a comparison of 
the colour-magnitude diagrams of the GCS of NGC 1399 in the Washington system 
and in V-I. The many unresolved background galaxies (a few stars are mixed in as 
well) with an excess in the blue filter C fall outside the colour-range of the globular 
clusters in C-Tl, whereas they populate that range in V-I and so add significantly to 
the background at faint magnitudes. 
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Fig. 3. This plot is based on HST observations of a sample of 14 early-type galaxies by 
Larsen et al. |56|. It shows the differences of the TOMs derived from the t5-function 
fits leaving the width either variable or non-variable versus the TOM derived by fits 
with a variable width. The two most deviating galaxies are NGC 1023 and NGC 3384. 
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least not in the case of well observed GCLFs (e.g. see Delia Valle et al. for 
the GCLF of NGC 1380). 

Moreover, the width of the adopted fitting function can be left free or can be 
fixed. Larsen et al. [HE] performed t5-function fits to their sample of 14 early-type 
galaxies both with the width as a fit parameter and with a fixed with of a t = 1.1 
mag. To gain an impression of the effect on the TOM, we show FigJ5J where the 
TOM corresponding to a free width is plotted versus the difference TOM(var)- 
TOM(nonvar). Leaving the two outliers NGC 1023 and NGC 3384 aside, the 
standard deviation of the differences is 0.13 mag. Then there are different ways 
to fit, for example maximum likelihood methods or direct fits. 

One must also not forget the uncertainty of the foreground absorption and 
another factor which is difficult to nail down: the photometric calibration of the 
respective data set and the actual realization of the used photometric standard 
system. 

The above error sources are always there, even if the TOM would be strictly 
universal, which one would not expect: For a given mass the luminosity of a GC 
depends on its metallicity in the sense that metal-poor clusters are brighter in 
the optical (e.g. Girardi |30|L but the metallicity distribution within a cluster 
system is not expected to be the same for all early-type galaxies. So one should 
principally correct for this as well (Ashman et al. 0). There is also empirical 
evidence for a metallicity dependence of the TOM: Larsen et al. jSE] find by 
HST observations in V and I of a sample of 15 early-type galaxies the TOMs 
for red clusters to be fainter than those for blue clusters by Amy « 0.4 mag, 
which is somewhat larger than predicted by theory. However, as we shall see, it 
is possible that part of this difference is due to an intrinsically fainter TOM of 
the red cluster population, so it is difficult to quantify the metallicity effect. 

Also if the initial cluster mass distribution would be the same in all galax- 
ies, destruction processes like disk shocking or evaporation are expected to act 
differently in different environments and may create intrinsically varying TOMs 
(see section 10 for this topic). Last, but not least, one has to assume that the 
members of a GCS all have the same old age, while we shall see that the number 
of examples where this is not the case is growing. 

Given all these possible error sources, it may come as a surprise that GCLFs 
seem to work so well as distance indicators, and it is plausible that an accuracy 
of, say, 0.2 mag or less can only be achieved in the case of rich GCSs and a high 
quality dataset. 

3 The Galaxy and M31 

The two massive galaxies where the GCLF can be best observed down to faint 
clusters are the Milky Way and M31. A calibration of the absolute TOM therefore 
via these galaxies always has to face the caveat that both are spiral galaxies 
and the application to early-type galaxies may not be justified. However, as we 
will see, the zero-point gained from using the Galaxy and M31 as fundamental 
calibrators is in very good agreement with the one obtained from the comparison 
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with the method of surface brightness fluctuations (Tonry et al. which at 

present offers the largest and most homogeneous catalog of distances to early- 
type galaxies. 

The history of the investigations of the galactic GCS and that of M31 involves 
the work of many people. To be short, we refer the reader to Harris and 
Barmby et al. 0] (and references therein) for the Galaxy and M31, respectively. 
Harris gSj quotes My = -7.40 ± O.f f for the galactic TOM and a = f .15 ± 
0.11. Barmby et al. 4 quote for the apparent TOM of the M31 system my = 
16.84 ±0.11 and a = 1.20 ±0.14. The distance modulus of M31 is m-M M31 = 
24.44 ± 0.2 (Freedman & Madore [23), which translates into an absolute TOM 
of My = -7.60±0.23. The weighted average of these two TOMs is -7.46±0.18, 
which we will compare with the distance moduli derived from surface brightness 
fluctuations. 

4 The data 

During the last few years, many new TOMs of early-type galaxies in the V-band 
have been published. The majority of them are based on HST observations and 
stem from the papers by Kundu & Whitmore [50|, Kundu & Whitmore |5*T] . 
and Larsen et al. JHE3- Other new papers on individual galaxies are from Okon 
& Harris [71] . Kavelaars et al. @7J, Woodworth & Harris |l()4j . Drenkhahn & 
Richtler ^H] • For TOMs published earlier we refer the reader to the compilation 
of Ferrarese et al. |20| and references therein. As mentioned above, the main 
problem with such a data set is its inhomogeneity for a variety of reasons. For 
example, Kundu & Whitmore [SO] and Kundu & Whitmore [HI] fitted Gaussians 
with both variable and fixed dispersions (1.3 mag) to their GCLFs. We adopt 
their TOMs resulting from the fixed dispersions because of the larger number 
of galaxies included, leaving out a few TOMs with very large uncertainties. 
Larsen et al. fitted t5-functions with both non- variable and variable widths, 
from which we adopt the latter because the scatter of the dispersions points to 
real differences. However, the TOMs are not strongly influenced by whether the 
dispersions are fitted or keep fixed. Since we are interested rather in the bulk 
properties of the available data than in hand-selected data according to certain 
quality criteria, we included also work which was mentioned but rejected by 
Ferrarese et al. 

We end up with 102 TOMs (corrected for foreground extinction and including 
a few double and triple measurements) in the V-band for 74 galaxies, which 
should be almost complete from the present day back to 1994. 

5 The Hubble diagram 

Can we say something about the Hubble constant from our data set, assuming 
that the TOM indeed has the universal value of My = —7.46 ± 18, adopted from 
the Milky Way and M31? A set of standard candles whose redshifts are only due 
to their recession velocities give a straight line in the Hubble diagram when their 
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log(c*z) 

Fig. 4. The upper panel shows all of our collected TOMs versus their recession veloc- 
ities, defined here as the radial velocities related to the microwave background. Many 
of these galaxies obviously have peculiar velocities of the same order as their recession 
velocities, in which case this definition is not adequate. The lower panel selects those 
galaxies with TOM uncertainties less than 0.3 mag and log(cz) larger than 3.2 to reveal 
a Hubble constant of 83 km/s/Mpc. Also here, the TOMs do not define very well a 
slope of 5 in the Hubble diagram. Note that the group at V=27 are not directly mea- 
sured TOMs, but deduced from surface brightness fluctuations, see Lauer et al. |6U| 
and Kavelaars et al. |47|. 

apparent magnitudes (their TOMs in our case) are plotted versus their redshifts 
according to 

m = 5 • log(c ■ z) - 5 • log{H ) + M - 25, 

where Hq is the Hubble constant in units of km/s/Mpc and M the constant 
absolute magnitude of the standard candles. 

The upper panel of Fig^ shows the Hubble diagram for our entire database. 
The velocities of the host galaxies have been individually related to the microwave 
background (which of course is not a good approach). It is obvious that such 



Globular Clusters 9 



a diagram is not suitable for deriving the Hubble constant. Many objects show 
radial velocities which are simply not in the Hubble flow, most strikingly for 
NGC 4406 (which is represented by the double measurement with the lowest 
velocity). If we select galaxies with log c ■ z > 3.2 and furthermore only those 
with uncertainties less than 0.3 mag, we end up with about 20 galaxies. If these 
galaxies are used to calculate the zero point in the Hubble relation, it gives a 
Hubble constant of 83 km/s/Mpc, adopting a TOM of My = —7.5 mag. It is 
clear that one cannot be content with this. Standard candles in a Hubble diagram 
should define a straight line with a slope of 5, whereas the slope in this diagram 
is clearly steeper. To resolve this discrepancy, one has to carefully look into each 
individual GCS, select those TOMs with the highest degree of trustworthiness, 
and then investigate the recession velocities of individual galaxies. The measured 
radial velocities of galaxies within the space volume under consideration may not 
be good indicators for their recession velocities due to the existence of large scale 
peculiar motions, which are under debate (e.g. Tonry et al. |9"T]1. 

Following Kavelaars et al. 0Zj , a better way might be to consider only groups 
of galaxies, average the TOMs and assign a recession velocity to each group. 
Kavelaars et al. use the Virgo, the Fornax and the Coma cluster and arrive at 
69 ± 9 km/s/Mpc for the Hubble constant. But to fix the recession velocities 
even for these three galaxy clusters is far from trivial. 

To avoid very lengthy discussions, a better way of deriving the Hubble con- 
stant is perhaps the use of standard candles which are so distant that peculiar 
velocities act only as minor perturbations of the Hubble flow, i.e. the Hubble 
diagram of Supernovae la (Freedman et al. |24p. 

6 The comparison with surface brightness fluctuations 

To evaluate the accuracy and reliability of the method of the GCLF, we must 
compare it with other distance indicators of early-type galaxies. This results 
in a complicated task, if one's objective is to select the most reliable measure- 
ments, to quantify possible biases inherent to different methods, and to discuss 
the uncertainties claimed by the authors. See for example Ferrarese et al. |2()| . 
who conclude that GCLFs do not provide reliable distances, mainly based on 
a deviating distance to the Fornax cluster, and Kundu & Whitmore J3U], who 
contrarily find GCLF distances as least as accurate as distances from surface 
brightness fluctuations. We do not want to follow these lines but rather investi- 
gate what can be seen from the entirety of TOMs if they are compared with a 
distance indicator which provides distances to most of our GCSs. 

Today, the most homogeneous and largest sample of distances to early-type 
galaxies is the catalog resulting from the survey of surface brightness fluctuations 
(SBFs) (Tonry et al. J§2|; see also the preceding papers by Tonry et al. [50] , 
Blakeslee et al. 0, and Tonry et al. which contains distances to about 300 
galaxies. Therefore we restrict ourselves to a comparison with this important 
distance indicator. Basically, it analyzes that part of the pixel-to-pixel scatter 
of a CCD image of an early-type galaxy which is caused by the finite number 
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Fig. 5. The upper panel shows the difference between the TOMs and the SBF distance 
moduli. For faint TOMs exists a trend to overestimate the distance with respect to the 
SBF distance. The lower panel selects those galaxies with uncertainties of both the 
TOM and the SBF distance less than 0.2 mag. The mean difference agrees very well 
with the absolute TOMs from the Milky Way and M31. 



of bright unresolved stars covered by each CCD pixel. These fluctuations of 
the surface brightness are large for nearby galaxies and small for more distant 
galaxies. 

Fig.|S]plots for all galaxies in our database (irrespective of whether there are 
double or triple measurements) the TOM versus its difference to the distance 
moduli from Tonry et al. • The error bars of the differences simply are the 
square roots of the quadratic sums of the uncertainties in the GCLF and SBF 
distance moduli. The first impression seems to be somewhat discouraging. Where 
we would have expected to see a horizontal line at an ordinate value of -7.5 with 
some scatter, we see a large spread with often dramatic devations, particular for 
the fainter TOMs. What is striking is that the deviating galaxies do not scatter 
symmetrically around a mean value, but that the faint TOMs give systematically 
larger distance moduli than do the SBFs. A direct and naive conclusion could be 
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that perhaps the very faint TOMs are observationally not reached and that an 
extrapolation from the bright end of the luminosity function to the TOM gives 
a TOM which is systematically too faint. In fact this is not the case and the 
strongly deviating TOMs belong to interesting galaxies (we come back to this 
point). 

But also at the bright end there are irritations. The deviating galaxy at -8.6 
is NGC 4565, and even the one with the brightest TOM, the Sombrero galaxy 
NGC 4594, does not fit very well to our assumed universal value. Both are the 
only spiral galaxies in our sample. We note that the GCS of NGC 4565 is very 
poorly populated (Fleming et al. |23), so this deviation might not bear much 
significance. 

7 Absolute TOMs and the distances to Virgo and Fornax 

However, if we select according to the quoted uncertainties, the situation starts 
to look better. The lower panel of Fig0 plots all galaxies where the uncertain- 
ties of both the SBF distance and the TOM according to the various authors 
are lower than 0.2 mag. The dispersion of the scatter is 0.25 mag and thus is 
compatible with the claimed selection. Thus we can confirm the statement by 
Kundu & Whitmore [S3 that the GCLF distances, at least for the sample under 
consideration, are not less accurate than the SBF distances. The mean differ- 
ence is -7.51 mag with a dispersion of 0.24 mag and thus in excellent agreement 
with the zero-points coming from the Milky Way and from M31. These three 
zero-points give a weighted mean of —7.48 ± 0.11. 

The average TOM of 8 galaxies in the Fornax cluster is 23.79 mag with a 
dispersion of 0.17 mag, the one for the Virgo cluster (16 galaxies) is 23.62 with 
a dispersion of 0.16 mag, which translate into distance moduli for Fornax and 
Virgo of 31.27 ± 0.2 and 31.10 ± 0.2, respectively. The corresponding distance 
moduli from the SBFs are 31.02 ± 0.15 and 31.49 ± 0.12. A discussion of the 
absolute calibration is not our objective. However, we can conclude that indeed 
many GCLFs can provide good distances but one is reluctant to label the GCLF 
"universal" at this point because there are too many deviations with the SBF 
distances. We shall see that these deviations are apparently related to the exis- 
tence of intermediate-age populations in early-type galaxies. 

8 Deviations and intermediate-age populations 

Like in human society, deviations from the norm may be sometimes more inter- 
esting and illustrative than reconciliation with it. Let's look at FigH3 Plotted 
are those TOMs which deviate from the "universal" TOM by more than what 
is suggested by their uncertainties. Since we see from Fig. [3] that the scatter 
is by no means symmetric but that the most striking deviations prefer a TOM, 
which is systematically fainter than expected from the SBF distances, only these 
fainter TOMs are shown. 
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Fig. 6. This graph shows those galaxies, whose GCLF distances deviate more than 
suggested by the measurement uncertainties from SBF distances. The upper panel 
shows the elliptical galaxies (t-parameter less than -4) , the lower panel the SO-galaxies 
(t-parameter higher than -4). In the majority of these objects, one finds evidence that 
the stellar populations are not exclusively old. A certain fraction of intermediate-age 
clusters would make the deviation from the SBF distance understandable, as explained 
in the text. 



Among the elliptical galaxies, the largest deviation is shown by NGC 3610, 
admittedly with a large error. But more interesting is the fact that this galaxy 
violates one important condition for the TOM to be a viable distance indicator, 
namely that its GC population is old. NGC 3610 is known to host GCs of 
intermediate-age. Whitmore et al. ||!|| and Whitmore et al. |1U2| estimate an age 
of about 4-6 Gyr for the red (and presumably metal-rich) GCs which plausibly 
have their origin in a merger event (Schweizer & Seitzer |80| ) . However, many of 
the metal- rich clusters probably had been brought in by the progenitor galaxies. 
Strader et al. 85 find among their sample of 6 metal-rich clusters only one with 
an age of 1-5 Gyr. But since the other clusters are located in the outer halo, this 
cannot strongly constrain the fraction of intermediate-age metal-rich clusters. 

As we later shall discuss, the GCLF is expected to change its TOM, resulting 
in fainter TOMs for younger cluster populations. Indeed, Whitmore et al. 102 
find a TOM of 25.44 ± 0.1 for the blue clusters alone, while no TOM is visible 
at all for the red cluster population. This decreases the difference to the SBF 
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distance, however, it still remains large. On the other hand, the existence of an 
intermediate-age population also influences the SBF distance by enhancing the 
fluctuation signal (mainly through the enhanced number of asymptotic giant 
branch stars) and thus would lead to a spuriously smaller distance without any 
correction term. Normally, the color V-I is used to correct for differences in the 
population (see Liu et al. and Blakeslee et al. [7] for a deeper discussion). 
Whether this correction is always sufficient or fails in some casesx, cannot be 
discussed here. 

So the suspicion arises that a difference between the GCLF distance and 
the SBF distance may in general be produced by intermediate-age GCs. This 
conjecture indeed gets support by looking at other galaxies in FigHj] Besides 
NGC 3610, younger clusters have been detected in NGC 4365 (Larsen et al. |59"] . 
Puzia et al. |7S], also conjectured to be a merger remnant (Surma & Bender |57)1. 
Note, however, that Davis et al. did not find evidence for intermediate-age 
populations from their integral-field spectroscopy in the galaxy iself. 

Other ellipticals show strong i^-lines, indicating as well a younger popula- 
tion, and NGC 4636 hosted a supernova la, whose progenitors should also be of 
intermediate age (e.g. Leibundgut |61|L 

However, there are also examples where the existence of an intermediate-age 
population is not supported by the present literature. For NGC 3379 and NGC 
4472, the difference to the SBF distance is anyway marginal, perhaps still so 
for NGC 4660. The case of NGC 4291 is hard to assess because of the large 
uncertainties, but NGC 4473 and NGC 5846 pose a problem. The spectroscopic 
evidence for intermediate-age populations normally come from the central re- 
gions and it may be that the outer parts, from which the globular clusters are 
sampled, still host younger populations. The other possibility is that either the 
SBF distance or the TOMs are erroneous. In any case one has to wait for further 
observations. 

Turning to the SO's, NGC 1023 and NGC 3384 are perhaps candidates for 
hosting intermediate-age populations. Both galaxies show indications of star for- 
mation activity in their inner regions (Kuntschner et al. [521, Sil'chenko (53]). 
However, the faint TOM of NGC 1023 does not seem to be related to intermediate- 
age clusters. Larsen & Brodie [55] identified beside the "normal" compact GCs 
(both red and blue) a population of faint extended red GCs. The inclusion of 
these latter objects in the luminosity function is mainly responsible for the de- 
viating position of NGC 1023. Leaving them aside results in a distance modulus 
well agreeing with the SBF distance. A similar finding is reported for NGC 3384 
by Larsen et al. [55] . Brodie & Larsen |S] found that these faint extended clusters 
belong to the disk populations of their host galaxies and quote an age of at least 
7Gyr. 

NGC 4550 contains two counterrotating stellar disks (Rix et al. (ZH)) and 
molecular gas has been detected by Wiklind & Henkel |l()3j which is supposed 
to have its origin in a recent accretion event. Similarly striking findings are not 
reported for NGC 1553 or NGC 1201. However, both are shell galaxies (e.g. 
Longhetti et al. HU), which hints at earlier interactions or mergers. 
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NGC 524 again is a supernova la host galaxy. The question whether the 
appearance of a supernova of type la in an early-type galaxy always indicates an 
intermediate-age stellar population is beyond the scope of this article, but a few 
remarks on GCSs of la host galaxies are appropriate. NGC 1316 in the Fornax 
cluster, a merger remnant and host to two SN la's, has a GCS where 2-3 Gyr old 
clusters have been found, probably formed during the merger event (Goudfrooij 
et al. EH], |2S]). Deep VLT and HST photometry does not reveal a TOM; the 
GCLF increases steadily down to beyond the observation limit (Grillmair et al. 
p?fT) , Gilmozzi, this conference) . In the work of Gomez & Richtler [2Sj who quote 
a TOM which is in good agreement with the SBF distance, the TOM was not 
actually reached, but extrapolated, and the agreement with the SBF distance 
perhaps stems from the fact that in the outer region, where this data has been 
sampled, the fraction of intermediate-age clusters is low. 

Other early-type la host galaxies with investigated GCSs, where intermediate- 
age cluster populations have been identified, are NGC 5018 (SN 2002 dj) (Hilker 
& Kissler-Patig gU) and NGC 6702 (SN 2002cs) (Georgakakis et al. [29 ). Un- 
fortunately, NGC 6702 is too far for an analysis of its GCLF and the TOM of 
the NGC 5018 system must be largely extrapolated, so it remains uncertain. 

But we also have examples of la hosts, where the GCLF distance agrees 
quite well or is even smaller than the SBF distance, e.g. NGC 4621 (2001A), 
NGC 1380 (1992 A), NGC 4526 (1994D) and NGC 3115 (1935B). If there are 
intermediate-age populations in these galaxies, they do not seem to contaminate 
the GCLF. 

An interesting note regarding la host galaxies and GCSs can be made from 
the paper of Gebhardt & Kissler-Patig [2U • These authors analyze the V-I colour 
distribution of the GCs of a sample of early- type galaxies. Their "skewness" 
parameter measures the asymmetry of the colour distribution with respect to 
the mean colour. The two GCSs which are skewed strongest towards red (e.g. 
metal-rich) clusters both belong to la host galaxies (NGC 4536, NGC 4374) as 
well as does the fourth in this sequence (NGC 3115). 

All this, of course, does not mean that in those cases where GCLF and SBF 
distances agree within the uncertainties, the stellar populations are necessarily 
old. However, a comparison with the compilation of galaxy ages by Terlevich & 
Forbes IHH] reveals that among the ellipticals, only NGC 720 (3.4 Gy) is quoted 
with an age lower than 5 Gyr. Among the SO's we have only NGC 3607 (3.6 
Gyr) and NGC 6703 (4.1 Gyr), i.e. strikingly less candidates for hosting younger 
populations than among the deviating ones. 

Summarizing, it seems that many of the cases where the SBF distance does 
not agree with the GCLF distance, can be related to the presence of intermediate- 
age populations, particularly among the ellipticals. 

Tab E] lists all galaxies in Fig. Elwith their TOM, its difference with the SBF 
distance, references for the TOM and a reference for other properties of the host 
galaxy. 
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Table 1. A list of galaxies whose TOMs indicated larger distances than the SBF 
distances beyond the uncertainty limits. For many of these galaxies one finds evidence 
for the existence of intermediate-age (IM) populations. 



Name 



TOM 



diff(TOM-SBF) Ref. 



Remarks 



Ellipticals 



N3610 
N5322 
N4589 
N0584 
N4636 
N4291 
N4697 
N5846 
N4458 
N4473 
N4660 
N4365 
N3379 
N4472 



26.49 
26.30 
25.22 
24.96 
24.10 
25.30 
23.50 
25.08 
24.20 
23.91 
23.39 
24.37 
22.82 
23.75 



±0.65 
±0.58 
±0.39 
±0.36 
±0.10 
±0.44 
±0.20 
±0.10 
±0.36 
±0.11 
±0.18 
±0.15 
±0.07 
±0.05 



-5.16 ±0.69 
-6.17 ±0.62 
-6.49 ± 0.45 
-6.56 ±0.41 
-6.73 ±0.16 
-6.79 ±0.54 
-6.85 ±0.24 
-6.90 ± 0.22 
-6.98 ±0.38 
-7.07 ±0.17 
-7.15 ±0.26 
-7.18 ±0.23 
-7.30 ±0.13 
-7.31 ±0.11 




IM clusters 
IM age 
Hj3 strong 
Hj3 strong 
la host 
old? 
IM age 
old? 

Hf3 strong 

old? 

old? 

IM clusters 

old? 

old? 



SO 



N3384 


23.30 ±0.13 


-7.02 ±0.19 


56 




extended GCs, Hf5 strong 


N4550 


24.08 ±0.16 


-6.92 ±0.26 






molecular gas, merger 


N0524 


25.00 ± 0.40 


-6.90 ± 0.45 






la host 


N1023 


23.53 ± 0.28 


-6.76 ±0.32 




MM 


extended GCs, IM nucleus 


N1201 


25.00 ± 0.60 


-6.53 ±0.67 




shell galaxy 


N1553 


25.20 ±0.60 


-6.14 ±0.62 


64 




shell galaxy 



9 Why does it work? 

What could be the reason for an universal TOM of old GCSs? A globular clus- 
ter with My = —7.5 mag has a mass of about 150000 Mq, adopting an average 
M/Ly of 2.5 (Pryor & Meylan [72] )• Is this particular mass somehow distin- 
guished? One has to realize that the magnitude scale is logarithmic. Binning 
in linear luminosity units instead of magnitudes, we would not see any striking 
feature at the luminosity corresponding to the TOM. After remarks by Surdin 
|5rj] , Racine [THj , and Richtler [7S| , regarding the power-law nature of the linear 
luminosity function of galactic GCs, McLaughlin [HS] put this concept on a for- 
mal basis. If the luminosity function can be described as NdL ~ L a ^ L \ where 
N is the number of clusters found in the luminosity interval L±dL and a(L) 
a function of L, then one has in mag nitudes NdM v ~ io°-«Mi-*) ; i e thc 
TOM is found where a(L) just changes from smaller than -1 to larger than -1. 
So the location of the TOM does not express a specific physical property at this 
particular mass. However, the underlying universal property must be a universal 
mass function. Harris & Pudritz 01] first investigated the mass function of GCSs 
of different galaxies, assuming a constant M/L. They found that such diverse 
systems as that of the Milky Way and of M87 can be described by a common 
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20 21 22 23 

mag(R) 

Fig. 7. This plot shows the luminosity lunction in the R-band for the GCS of NGC 
1399, the central galaxy in the Fornax cluster (Dirsch et al. |13|). This luminosity 
function comprises about 2600 clusters brighter than R=23 and thus is one of the best 
available. For magnitudes fainter than R=20.5, the linear luminosity function is well 
represented by a power-law with an exponent of about -2 (the slope s in the present 
diagram is related to the power-law exponent a by a — s/0.4 - 1). It steepens consid- 
erably for brighter magnitudes. However, a representation by a log-normal function for 
the entire magnitude range is as good. 

power-law exponent of a « —1.8 for masses higher than about 10 5 solar masses. 
Larsen et al. •">(> found in their larger sample on the average a = —1.74 ± 0.04 
between 10 5 and 10 6 solar masses. However, in very rich GCSs, such as that of 
M87 or NGC 1399, the slope becomes distinctly steeper for cluster masses larger 
than about 10 6 solar masses (see FigJTJ. 

Ten years ago, GCSs had been almost exclusively associated with old stellar 
populations. Meanwhile, systems of young globular clusters have been detected 
in many merging galaxies, the most prominent ones being the Antennae NGC 
4038/4039 (Whitmore & Schweizer .98 ) and NGC 7252 (Whitmore et al. 1993) 
(see Whitmore |101) for a complete listing until 2000), but also in normal spiral 
galaxies (Larsen & Richtler Larsen & Richtler |54j ) . 

Determinations of the luminosity functions resulted so far consistently in 
power-laws with an exponent of about -2, without compelling indications that 
this exponent changes over the observed luminosity range as in the case of the 
GCSs of giant ellipticals (Whitmore |101) 1. given the uncertainties caused by 
internal extinction and by the age spread among a cluster system. There was 
some debate regarding the mass function of GCs in the the Antennae as de- 
rived from the luminosity function. The Antennae may show a bend at about 
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My = —11, becoming steeper towards the bright end (Whitmore |101J . Zhang 
& Fall |105p . Fritze-v. Alvensleben [251 ) found a log- normal mass distribution 
like for old systems, which was contradicted by Zhang and Fall ( |1()5| ). who at- 
tributed this difference to the effect of varying extinction and ages, and found a 
uniform power-law. If we assume that young GCs are born obeying a universal 
luminosity function like dN/dL ~ L~ 2 , and accordingly with a mass function 
of the same shape, then we must ask, what processes can transform such a lu- 
minosity function into the approximately log-normal luminosity functions of old 
GCSs. If these processes work in a universal manner, then the universality of 
the TOM could be explained. 

10 How does an initial cluster mass function change with 
time? 

3 clusters with 1,2,4 10**5 solar masses 
R a = 6 kpc, R_p = 4 kpc 




only relaxation 

relaxation 
+ grav. shocks 

relaxation 
+ grav. shocks 
+ stellar evol. 



t (Gyr) 

Fig. 8. This plot is taken from Fall & Zhang |19| . It shows for three different clusters 
(10 5 ,2 • 10 5 ,4 ■ 10 5 solar masses) the cumulative growth of the relative mass loss by 
two-body relaxation (upper line), two-body relaxation plus tidal shocks (middle line), 
and added to that the mass loss from stellar evolution (lower line). 

A young star cluster is exposed to different destruction mechanisms. If it is 
still young, mass loss from massive star evolution plays an important role. At 
later times, two-body relaxation, dynamical friction, and tidal shocks, when the 
cluster enters the bulge region of its host galaxy or moves through a disk, can 
be efficient in decreasing the cluster's mass, depending on its mass, it's density 
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and it's orbit in the host galaxy. The most general statement is that low- mass 
clusters are more affected by disruption processes than high-mass clusters, so 
an initial power-law of the mass distribution is more strongly distroyed on the 
low mass end and may develop a shape which finally resembles a log-normal 
distribution. 




log (M/M @ ) 

Fig. 9. This plot is taken from Fall & Zhang |19|. It shows the evolution of a cluster 
system for two different initial mass functions: a power-law with an exponent of -2 
(upper panel) and a Schechter function (lower panel). 



Many people have worked on this problem, among them Aguilar et al. pQ, 
Okazaki & Tosa [ZD], Elmegreen k Efremov [nj, Gnedin & Ostriker [H], Murali 
& Weinberg (1997a,b,c) Vesperini [M|,[55], Fall & Zhang [TU]. We cannot present 
all work in detail, instead we choose the analytical model by Fall & Zhang in 
order to illustrate the most important results. Fig. [5] (Fig. 1 of Fall & Zhang 
2001) shows the time evolution for three different masses for a cluster which is 
on a slightly elongated orbit. The dotted lines indicate the effect of two-body 
relaxation alone. The dashed lines additionally include gravitational shocks, and 
the solid lines add the effect of mass loss by stellar evolution. 
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Under a wide variety of conditions, the mass function of a GCS develops 
a peak which is progressively shifted to higher masses, as the evolution of the 
cluster system proceeds. After, say, 12 Gyr, Fall & Zhang get from their model a 
peak mass (in logarithmic bins) which may well represent the mass corresponding 
to the TOM observed in the Milky Way or in elliptical galaxies (FigEJ. 

However, it seems that the assumption of a power-law with an exponent 
around -2, as suggested by the young cluster systems in merging galaxies, cannot 
reproduce well the log-normal shape in the mass-rich domain observed in many 
galaxies. This is because the shape of the mass function above a few times 
10 6 Mq practically does not change by evolutionary processes. Instead, an initial 
log-normal mass function works much better in resembling the bright end of the 
luminosity function of ellipticals (Vesperini , [96 ) (but see the section on the 
brightest clusters). 

The dynamical evolution of a GCS may raise doubts on the general quality of 
the GCLF as a distance indicator, if the evolutionary history of a GCS is not neg- 
ligible. The GCLF might also depend on whether the TOM is measured at small 
or large galactocentric radii. In the inner parts of a galaxy, the TOM is expected 
to be brighter. Gnedin jjH] finds significant differences in this sense for the Milky 
Way, M31, M87, which for M31 has been confirmed by the improved sample of 
Barmby et al. [I]. Also the brighter TOM (with respect to the SBF distance) of 
the Sombrero may have its explanation in the dynamical history of this GCS. 
The Sombrero possesses an extraordinary large bulge, where dynamical shocks 
might work more efficient than in other galaxies (naively assuming, of course, 
that the SBF distance is correct). Note, however, that the HST-observations by 
Larsen et al. U2] reveal a GCLF for the Sombrero whose TOM is not very well 
defined. 

One of the best investigated galaxies among those which shows a marked 
difference between the GCLF and the SBF distance is the elliptical galaxy NGC 
3610. Scorza & Bender |^ found a disk and other morphological signatures 
indicating previous interaction or a merger event. It's location in Fig0 corre- 
sponds to the TOM quoted by Kundu & Whitmore [SOj- In a subsequent paper, 
Whitmore et al. |l()2j performed a more detailed investigation of the GCS of 
NGC 3610, based on new HST data. Fig. I1UI shows the LFs separately for the 
blue and the red clusters. While for the blue clusters the TOM is measured to be 
at V = 25.44±0.1, the red clusters show a LF rising until the photometry limit. 
Whitmore et al. combined the destruction model of Fall & Zhang with evolution- 
ary models of stellar populations. The resulting model LFs are indicated in the 
lower panel. The data are not yet deep enough to show a turn-over for the red 
clusters, which by the models is predicted to be at around V ~ 26. Whitmore et 
al. state that in the context of the Fall & Zhang models, the brightening of the 
TOM during the dynamical evolution of the cluster system is almost completely 
balanced by the fading of the stellar population during this time. This may well 
be an explanation for the universality of the TOM. However, given the approx- 
imate nature of the analytic models of Fall & Zhang and the the dependence 
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of the destruction processes on the actual environment, this probably does not 
apply to every galaxy. 

We therefore can conclude that in the case of NGC 3610, a large part of 
the deviation of the GCLF distance from the SBF distance comes from the fact 
that the contribution of the presumably younger red clusters causes a fainter 
TOM than from the blue, metal-poor and presumably older clusters alone. But 
also when we use only the TOM of the blue clusters to determine the distance 
modulus, which then would be 32.94, a significant difference remains to the SBF 
modulus, which is 31.65. This cannot be resolved here. The modelling of SBF's 
accounts for the population structure (Liu et al. Blakeslee et al. [Z]) but 
may fail in extreme cases. 

11 The brightest clusters 

Regarding the significance of the GCLF as a distance indicator, its shape in the 
domain of the brightest clusters is less important. But since dynamical models 
indicate that the GCLF for clusters more massive than about 10 Mq is not 
modified by destruction processes, they bear potential information about the 
formation of a GCS. Two different views on a GCLF like that of Fig0 exist: It 
can be seen as a power-law with an exponent around -2 with a cut-off at higher 
masses or it can be seen as a log-normal function. 

Adopting the first view, the GCLF would resemble in large parts the LF found 
in young cluster systems. The cut-off at high masses may have different reasons. 
Possibly these very massive clusters (the brightest clusters in NGC 1399 have 
about 10 7 M Q ) are not globular clusters in the normal sense, but the dynamically 
stripped nuclei of dwarf galaxies. In the Milky Way, we have uj Centauri as a 
possible example (e.g. Hilker & Richtler @S|)- I n this case, the LF at the bright 
end would depend on the accretion rate of dwarf galaxies, which pausibly is 
highest in massive galaxies in dense environments like NGC 1399 or M87. Or 
the formation history of the most massive clusters is principally different from 
the less massive ones. The peculiar cluster in NGC 6946 (Larsen et al. |5Z|> |5B|) 
with a mass of about 1O 6 M and an age of 15 Myr is surrounded by a round, 
star forming complex of about 600 pc diameter, which gives the impression of a 
disk-like structure with the massive cluster near its center. Such a configuration 
suggests that the cluster mass is determined or partly determined by accretion 
from a larger region, resulting in a steeper mass function for massive clusters. 

Taking the second view of a log-normal function, one has the possibility to 
relate such a shape to coagulation processes by which GCs might have been 
formed through the merging of smaller subunits. Based on ideas by Harris & 
Pudritz 01] , McLaughlin & Pudritz [HfJ developed a model, in which GCs form 
inside the cores of supergiant molecular clouds. These cores are built up by 
internal collisions and subsequent coagulation of smaller clouds. Star formation 
tends to partly disrupt these cores and in an equilibrium between coagulation 
and disruption, a mass spectrum of cores results, which directly resembles the 
GC mass spectrum. This is because the formation of GCs in these cores must 
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Fig. 10. This plot is taken from Whitmore et al. 102 . It shows the luminosity functions 
(LFs) of blue (upper panel) and red (lower panel) globular clusters in NGC 3610. The 
blue clusters exhibit a TOM at V = 24.55. The solid line marks a Gaussian fit with 
a = 0.66, the dotted line with a = 1.1, which obviously do not fit. The LF of the 
red clusters increases with no sign of a flattening. The thick solid line is a power-law 
fit with a = —1.78. The other curves are model LFs based on Fall & Zhang |19) in 
combination with population synthesis models of Bruzual & Chariot (unpublished). 
The thin solid line is the zero-age LF (power-law with a = —2), the others correspond 
to ages of 1.5 Gyr, 3 Gyr, 6 Gyr, and 12 Gyr (from top to bottom). The TOMs 
are hardly distinguishable because in these models the brightening of the TOM by 
dynamical evolution is balanced by the fading due to stellar evolution. 

occur with a high star formation efficiency in order for the GC to stay bound, 
i.e. the mass of the core is closely related to the mass of the final cluster. 

See also Burkert & Smith [S] who argue that the mass spectra of GCSs can 
be fitted with a form dN /dm ~ m~ 3 / 2 • exp(—m/m c ), where m c is a "truncation 
mass". Such a shape resembles the long-time solution of the coagulation model 
of Silk & Takahashi [SJ , initially starting with small progenitor clouds of equal 
mass. 

Although the relation between the GC mass spectrum and the mass spectrum 
of the progenitor clouds is open to speculation, the power-law interpretation of 
the GCLF has some attractive features over the log-normal law interpretation.lt 
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relates the GCLF of old clusters system with young ones, and it offers a simple 
explanation by a direct link to the mass spectrum of molecular clouds. 

It is amazing that molecular clouds in the Galaxy exhibit a mass spectrum 
resembling so closely that of GCs. See the introductory part of Elmegreen [T%| 
for a compilation of references. This power-law behaviour may be the result of a 
fractal structure of the interstellar gas caused by turbulence and selfgravitation 
(Fleck [22 j Elmegreen & Falgarone |17|. Elmegreen JHj). Therefore the univer- 
sality of the GCLF probably has its ultimate explanation in the universality of 
the interstellar gas structure. 

12 Conclusions 

We have seen that the method of globular cluster luminosity functions (GCLFs) 
allows one to determine distances to early-type galaxies, which are as accurate 
as those derived from surface brightness fluctuations (SBFs), once the condi- 
tions of high data quality and sufficiently rich cluster systems are fullhlled. The 
achievable accuracy of distance moduli is of the order 0.2 mag. The absolute 
turn-over magnitudes (TOMs), if calibrated by SBF distances, agree very well 
with those of the globular cluster systems of the Milky Way and the Andromeda 
nebula. Therefore the TOM is indeed a universal property of old globular cluster 
systems. 

The comparison of SBF distances with GCLF distances reveals however many 
discrepant cases, in which the GCLF distances are systematically larger than the 
SBF distances beyond the limits given by the uncertainties. In some elliptical 
galaxies, direct evidence for the existence of intermediate-age globular clusters 
is available. In others, intermediate-age stellar populations are indicated by a 
variety of findings, which again may suggest a certain fraction of intermediate- 
age globular clusters as well. The SO-galaxies NGC 1023 and NGC 3384 exhibit 
a population of faint extended red clusters, which cause a fainter TOM, if they 
arc included in the luminosity function. 

That a globular cluster system, consisting mainly of old clusters, in which 
intermediate-age clusters are mixed in, exhibits a fainter TOM, can be under- 
stood by the dynamical evolution of cluster systems. Young globular clusters, 
which are found in large numbers in merging galaxies, are formed according to a 
power-law mass function with an exponent around -2. The cluster system then 
undergoes a dynamical evolution where the mass loss of individual clusters is 
caused by two-body relaxation, tidal shocks and mass loss by stellar evolution. 
This results in a preferential destruction of low-mass clusters, modifying the 
initial power-law mass function in such a way that the corresponding luminos- 
ity function on the magnitude scale shows a TOM which becomes brighter as 
the system evolves. The fading of clusters by stellar evolution counteracts this 
brightening to some degree. 

The universality of the GCLF probably has its origin in the fractal struc- 
ture of the interstellar medium, which results in a power-law mass spectrum for 
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molecular clouds with an exponent of -2, similar to that found for young globular 
cluster systems. 
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